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In this lecture we discuss correlations of the QCD Dirac eigenvalues. We find that 
below a scale of E c ~ A/L 2 they are given by chiral Random Matrix Theory. This 
follows from analytical arguments based on partially quenched Chiral Perturbation 
Theory and is substantiated by lattice QCD and instanton liquid simulations. 



1 INTRODUCTION 

The low-energy limit of QCD is characterized by two parameters, the chiral 
condensate, E, and the pion decay constant, F. The chiral condensate is di- 
rectly related to the density of the smallest eigenvalues of the Dirac operator 
by means of the Banks-Casher formulae!. At finite volume, the pion decay con- 
stant determines the quark mass scale below which the QCD partition function 
is given by the zero momentum component of the chiral LagrangianBfcl. Be- 
low this scale the Dirac eigenvalues are constraint by Leutwyler-Smilga sum 
rules a On the other hand, the accumulation_pf small eigenvalues for broken 
chiral symmetry requires strong interactions □. From the study of complex 
systems (see Guhr et al. □ for a comprehensive review) we expect that mi- 
croscopic correlations of eigenvalues of such disordered system are universal 
and can be described by the simplest model in its universality class. Such a 
model is chiral Random Matrix Theory (chRMT). The interpretation is that 
the classical motion of quarks in Euclidean QCD with an additional artificial 
time dimension is chaotic. 

These ideas can be formulated more precisely in terms of partially quenched 
Chiral Perturbation TheoryO (pqChPT), which, in addition to the usual quarks, 
contains valence quarks and their superpartnersQ. If both the Vafa-Witten the- 
orem □ and the Goldstone theorem apply to this theory, only one mass term 
consistent with the chiral symmetries of the partition function can be written 
down. On the other hand, it can be shown □ that the zero-momentum sec- 
tor of pqChPT reproduces identically the chRMT distribution of the smallest 
eigenvalues up to a scale which is, the equivalent of the Thouless energy in 
the theory of mesoscopic systems Ejlj. In units of the level spacing, this scale 
is given by_p c = F L /n for a box of length L. This argument proves the 
conjecture t3 that the correlations of the smallest QCD Dirac eigenvalues are 
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our confidence in this_r.esult. These arguments have 
by instanton liquidEL!rE3 and lattice QCD simulations 

Because of the Banks-Casher relation, E = Trp(Q)/V, the eigenvalues near 
zero are spaced as l/p(0) = ir/EV. In order to study the approach to the 
thermodynamic limit it is natural to introduce the microscopic limit injsihich 
u = WE is kept fixed for V — > oo and the microscopic spectral densityEj 
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Our claim is that for u <C n c the microscopic spectral density is given by 
chRMT. Indeed, the chRMT result for (g) can be derived from pqChPTEl 

The chRMT partition function is introduced in section 2. Its domain of 
validity is discussed in section 3 and universality arguments are presented in 
section 4. Lattice results are shown in section 5. Before concluding we finish 
with applications of chRMT to QCD at nonzero chemical potential (section 6). 



2 CHIRAL RANDOM MATRIX THEORY 

In this section we will introduce an instanton liquid@@ inspired chiral RMT 
with the global symmetries of the QCD partition function but otherwise Gaus- 
sian random matrix elements. For JVjfiflavors m the sector with topological 
charge v such chRMT is defined byE30 

N f 

Zf ffiU (m 1 ,...,m Nf )= / DWl[det(V + m f )e-^™ tw , (2) 
J f= i 

where 

and W is a n x m matrix with v = \n—m\ and N = n+m. The matrix elements 
of W are either real ((3=1, chiral Gaussian Orthogonal Ensemble (chGOE)), 
complex ((3 = 2, chiral Gaussian Unitary Ensemble (chGUE)), or quaternion 
real ((3 = 4, chiral Gaussian Symplectic Ensemble (chGSE)). As is the case in 
QCD, we assume that v does not exceed y/N, so that, to a good approximation, 
n = N/2. The value of the Dyson index f3 is determined by the anti- unitary 
symmetries of the Dirac operator. If there are no anti-unitary symmetries 
the value of f3 = 2. If the square of the anti- unitary symmetry operator is 
1 we have /? = 1 and if its square is —1 we have (3 = 4. Together with the 
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Dyson ensembles, the GOE, the GUE and the GSE, the chiral ensembles can 
be classified in terms symmetric spaces cB. 

In this model chiral symmetry is broken spontaneously with chiral conden- 
sate given by S = limAr^oo 7rp(0)/N, where N is interpreted as the (dimen- 
sionless) volume of space time. For complex matrix elements (j3 = 2), which 
is appropriate for QCD with three-or more colors and fundamental fermions, 
the symmetry breaking pattern isL3 SU(Nf) x SU(Nf)/ SU(Nf). The average 
spectral density that can be derived from (Q) has the familiar semicircular 
shape. The microscopic spectral density for the chGUE is given byEIrEJ 

PS(Z) = I (J 2 a (z) ~ Ja+l(z)Ja-l(z)) , (4) 

where a = Nf + \ v\. The spectral correlations in the h ulk of the spectrum are 
given by the invariant random matrix ensemblesLj'Ej. 



3 DOMAIN OF VALIDITY OF chRMT 

The domain of validity is best discussed within the context of partially quenched 
Chiral Perturbation Theory (pqChPT). This is an effective field theory for the 
low-energy limit of a QCD like theory that in addition to Nf sea quarks con- 
tains valence quarks and their superpartners. It allows us to calculate the 
valence quark mass dependence of the chiral condensate defined as 

1 C 2m 
^ {mv) = -j dXm) —^-. ( 5 ) 

The spectral density follows from the discontinuity of S t ,(m tJ ), 

y{p(X)) = lim E(iA + e) - E(iA -e) = y ]>>( A + A*)), (6) 

k 

where the average (• • •) is with respect to the distribution of the eigenvalues. 
Similarly, the two-point spectral correlation function follows from the double 
discontinuity of the scalar susceptibility. 

In both cases we can identify an important scale where the inverse mass 
of the Goldstone modes corresponding to the valence quark is equal to the size 
of the box. Using the relation M = (m + m')S/F 2 , where F is the pion r 
constant, we find from ML = 1 that in terms of m this scale is given byt 

F 2 

E ° - Elf ^ 
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For m v -C E c we have shown that the valence quark mass dependence is given 
by the chRMT. The asymptotic result for the chRMT two-point correlation 
function follows from the double Goldstone pole in the scalar susceptibility. 
The conclusion is that if pqChPT describes correctly the low energy limit of 
QCD we have shown that the correlations of the Dirac eigenvalues close to zero 
are given by chRMT. 

Such picture is well-known from mesoscopic physics i. In this context E c 
is defined as the inverse tunneling time of an electron through the sample and 
is given by E c = hD/L 2 , where D is the diffusion constant. Another scale that 
enters in these systems is the elastic scattering time r e . Based on these scales 
one can distinguish three different regimesllil for the energy difference, 5E, that 
enters in the two-point correlation function: i) the ergodic regime for 5E <§C E c , 
ii) the diffusive domain for E c <C SE <C h/r e and iii) the ballistic regime for 
Se h/r e (not discussed below). For time scales corresponding to the ergodic 
regime an initially localized wave packet covers all of phase space. In this 
domain the eigenvalue correlations are given by RMT. In the diffusive domain 
an initially localized wave packet explores only part of the phase space resulting 
in weaker correlations between the eigenvalues. For earlier, .applications of 
localization theory to the chiral phase transition we refer tocireHl. 

Based on these ideas we can interpret the Euclidean Dirac spectrum as 
the energy levels of a Hamiltonian conjugate to_an additional artificial time 
dimension. According to the Bohigas conjecture c3 the eigenvalue correlations 
are given by RMT if and only if the corresponding classical motion is chaotic. 
We thus conclude that the classical time evolution of quarks in the Yang-Mills 
gauge fields is chaotic. 
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Figure 1: The number variance E2(n) versus n measured from an interval starting at A = 0. 
The total number of instantons is denoted by TV. 
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These ideas can be tested by means of lattice QCDBE3 and instanton liquidH 
simulations. In Fig. 1, we show E3 T*2[ri), defined as the variance of the 
number of levels in an interval containing n level on average, versus n for 
eigenvalues obtained from the Dirac operator in the background of instanton 
liquid gauge field configurations. The chRMT result, given by the solid curve, is 
reproduced up to about two level spacings. In units of the average level spacing, 
A = l/p(0) = n/T,V, the energy E c is given by n c = E c /A = F 2 L 2 /ir. For an 
instanton liquid with instanton density N/V = 1 we find that n c w 0.07y/~N. 
We conclude that chRMT appears to describe the eigenvalue correlations up 
to the predicted scale. 

The additional two flavors in Damgaard's relation E3 between the micro- 
scopic spectral density and the finite volume partition function arise naturally 
in pqChPTH, namely one valence flavor and its supersymmetric partner. 



4 UNIVERSALITY IN CHIRAL RANDOM MATRIX THEORY 

The aim of universality studies is to identify observables that are stable against 
deformations of the random matrix ensemble. Not all observables have the 
same degree of universality. For example, a semicircular average spectral den- 
sity is found for random matrix ensembles with independently distributed ma- 
trix elements with a finite variance. However, this spectral shape does not 
occur in nature, and it is thus not surprising that it is only found in a rather 
narrow class of random matrix ensembles. What is surprising is that the 
microscopic spectral density and the microscopic spectral correlators are sta- 
ble with respect to a much larger class of deformations. Two different types of 
deformations have been considered, those that maintain the unitary invariance 
of the partition functions and those that break the unitary invariance. 

In the first class, the Gaussian probability distribution is replaced by 
P{W) ~ exp(-A^fci a fcTr(W^ffi) fc ). For a potential with only a x and a 2 
different from zero it was shown L3 that the microscopic spectral density is 
independent of a 2 . A general proof valid for arbitrary potential and all cor- 
relation functions was given by Akemann et al. lii The essence of the proof 
is a remarkable generalization of the identity for the Laguerre polynomials, 
linin^oo L n (x/n) — Ja(2yfx) , to orthogonal polynomials determined by an 
arbitrary potential. It was proved from the continuum limit of the recursion 
relation for orthogonal polynomials. 

In the second class, an arbritrary fixed matrix is added to W in the Dirac 
operator (^) . It has been shown that the microscopic spectral density and the 
microscopic spectral correlations remain unaffected l 14 H 15 l for parameter values 
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that completely modify the average spectral density. 

Microscopic universality for deformations that affect the macroscopic spec- 
tral density implies the existence of a scale beyond which universality-breaks 
down. It can be interpreted naturally in terms of the spreading widtha. 

Based on the general form of of the pqChPT partition function one could 
argue that universality follows from its relation with chRMT. However, uni- 
versality can be formulated as the stability of the effective partition function 
with respect to variations of the distribution of matrix elements. The stability 
of the saddle-point manifold was demonstrated for the Unitary ensembles LJ. 



5 LATTICE QCD RESULTS 

In this section we consider correlations of lattice QCD Dirac eigenvalues. For 
a Dirac operator with a £^(1) symmetry the eigenvalues occur in pairs ±A&. 
Therefore we have to distinguish two different regions: the region near zero 
virtuality and the bulk of the spectrum. The Ua(1) symmetry is absent for 
the Hermitean Wilson Dirac operator which includes an additional 75 matrix. 

The relevant class of random matrix ensembles is determined by the anti- 
unitary symmetries of the Dirac operator. For the SU(2) color group, the 
anti-unitary symmetri£s_of the Kogut-Susskind (KS) and the Wilson Dirac 
operator are given byC3c3, 

[D ks ,t 2 K} = 0, and [ lh D w ^CKt*] = 0. (8) 

Because (T2K) 2 = —1 and (d^CKr-i) 2 — 1, we have (3 = 1 for Wilson fermions 
and (3 = 4 for KS fermions. This results in the chGUE for KS fermions and, in 
absence of the Ua{\) symmetry, in the GOE for Wilson fermions c3. For three 
or more colors there are no anti-unitary symmetries and the relevant ensembles 
are the chGUE and GUE for KS and Wilson fermions, respectively. 

The eigenvalue correlations, which are separated from the average spectral 
density by calculating them for the unfolded eigenvalues, are studied by means 
of moments of the number of levels, n^, in an interval containing n levels on 
average. Below we discuss the variance, ^(n) and the first two cumulants 
71 (n) and 72(H). The moments can be obtained either by ensemble averaging 
or by spectral averaging. Under the assumption of spectral ergodixjity, both 
procedures should give the same results. However, recent resultsEj indicate 
that the domain of RMT correlations is the much longer for spectral averaging 
than for ensemble averaging. This interesting observation certainly deserves 
more attention. 

Results obtained by spectral averaging for £2(71), 71 {n) and 72(71) for both 
KS and Wilson fermions with N c = 2 (see Fig. S) show an impressive agree- 
ment with the RMT predictions. This agreement extends to a wide range of 
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/3 values ranging from strong^pupling to weak coupling. The simulations for 
KS fermions were performed lJ for 4 dynamical flavors with ma = 0.05 on a 
12 4 lattice. The simulations for Wilson fermions were done for two dynamical 
flavors on a 8 3 x 12 lattice. Recently, impressive results were obtained fox. the 
nearest neighbor spacing distribution of N c = 3 staggered Dirac spectra Ea. 

Spectral ergodicity cannot be exploited in the study of the microscopic 
spectral density. In order to gather sufficient statistics, a large number of 
independent spectra is required. One way to proceed is to use instanton-liquid 
configurations which can be generated cheaply. In this case we find El that the 
microscopic spectral density is given by the chGOE and the chGUE for N c = 2 
and iV c = 3, respectively. Recently, this analysis was performed for lattice QCD 
Dirac eigenvalues!^. Results for 1416 quenched SU (2) Kogut-Susskind Dirac 
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spectra on a 10 4 lattice are shown in Fig. [| We show both the distribution of 
the smallest eigenvalue (left) and the microscopic spectral density (right). The 
results for the chGSE are represented by the dashed curves. We emphasize that 
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Figure 3: The distribution of the smallest eigenvalue (left) and the microscopic spectral 
density (right) for two colors and = 2.0. 

the theoretical curves have been obtained without any fitting of parameters. 
The input parameter, the chiral condensate, is derived from the same lattice 
calculations. The above simulations were performed at a relatively strong 
coupling of (3 — 2, but agreement with the-daGSE predictions was also found 
for (3 = 2.2 and for (3 = 2.5 on a 16 4 latticeEl 

In the case of two fundamental colors the continuum theory and Wilson 
fcrmions are in the same universality class. It is an interesting question of how 
spectral correlations of KS fermions evolve in the approach to the continuum 
limit. Certainly, the Kramers degeneracy of the eigenvalues remains. However, 
since Kogut-Susskind fermions represent 4 degenerate flavors in the continuum 
limit, the Dirac eigenvalues should obtain an additional two-fold degeneracy. 
We are looking forward to more work in this direction. 



6 APPLICATIONS OF chRMT AT /i ^ 

In the continuum formulation of QCD the chemical potential enters in the 
QCD partition function by the addition of the term to the anti-Hermitean 
Dirac operator, i.e. T> — > T> + /i7o- In a suitable chiral basis in which the 
matrix elements of (0jj|To|^iL) = $ki, the modification in the random matrix 
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partition function (g) corresponds to replacing the Dirac matrix T> by 

V=( + ° iW + ») (9) 

The term /ijo has the same anti- unitary symmetries as 7o<9o and therefore does 
not alter the correspondingjdassification of the Dirac operators. 

As is the case in QCDca, a nonzero chemical potential violates the anti- 
Hermiticity of the Dirac operator and its eigenvalues are scattered in the com- 
plex plane. One important result obtained from the chRMT at fi ^ is that 
the quenched approximation is the limit Nf — &_Q of the partition function with 
the determinant replaced by its absolute valueej. By replacing Li+iT one 
obtains a Landau-Ginzburg functional for the order parameter l3. A .similar 
Landau-Ginzburg functional has been derived from the NarnhiLjpDdelEiJ. For 
fx = this model undergoes a second oxder phase transition E3o'Ej. At T = 
we find a first order phase transitioned. A tricritical point is found for at 
m = 0, /i.^ and T ^ 0. For more discussion of this model we refer to the 
literatureEl 

2.0 — 

y I 
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Figure 4: Scatter plot of the real (x), and the imaginary parts (y) of the eigenvalues of the 
random matrix Dirac operator at fi ^ 0. 

We close this lecture with identifying possible universal behavior at /i ^ 0. 
Quenched numerical simulations have been performed for all three classes of 
(3. A cut along the imaginary axis below cloud of eigenvalues was found in 
instanton liquid simulations l£| for N c =_2 at fj, ^ (f3 = 1). In lattice 
QCD simulations with staggered fermionsEJ for N c = 2 (/3 = 4) a depletion of 
eigenvalues along the imaginary axis was observed whereas for N c = 3 (/? = 2), 
the eigenvalue distribution did not show any pronounced features. 
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In the quenched approximation, the spectral properties of the random ma- 
trix ensemble can be studied numerically by simply diagonalizing^a set of 
matrices with probability distribution (^). In Fig. |] we show results E3 for the 
eigenvalues of a few 100 x 100 matrices for \i = 0.15 (dots) of the quenched 
random matrix ensemble (^|). The solid curve represents the analytical result 
for the boundary of the domain of eigenvalues oE3. We observe the accumu- 
lation and depletion that was found in the previously mentioned simulations. 
This depletion can be understood as follows. For fi = all eigenvalues are 
doubly degenerate. This degeneracy is broken at /i =/= which produces the 
observed repulsion of the eigenvalues. 

The number of purely imaginary eigenvalues for (3—1 appears to scale 
as \fN which has been understood analytically for ensembles without a chiral 
structure E3. Obviously, more work has to be_done in order to arrive at a 
complete characterization of universal featuresEj of nonhermitean matrices. 

7 CONCLUSIONS 

We have presented both analytical and numerical arguments showing that the 
correlation of the QCD Dirac eigenvalues below a scale of ~ Kqcd/L 2 are 
given by chRMT. We conclude that if partially quenched chiral perturbation 
theory is the correct description of the low energy limit of QCD, we have 
shown that below the Thouless energy, the eigenvalue correlations are given by 
chRMT. Our confidence in these arguments has been strengthened by univer- 
sality studies. Universal features in Dirac spectra at nonzero chemical potential 
have been identified. 
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